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Abstract;  The  solubility  of  CL20  in  supercritical  carbon  dioxide  (CO2) 
was  evaluated  using  the  Peng-Robinson  cubic  equation  of  state.  Thermo¬ 
dynamic  properties  of  CL20  were  required  for  this  estimation.  A  compre¬ 
hensive  search  of  the  literature  revealed  little  data  regarding  the  critical 
point  of  the  compound.  Also  lacking  were  vapor  pressure  data  for  the 
compound.  Critical  properties  as  well  as  other  required  thermodynamic 
properties  were  estimated  using  a  variety  of  available  estimation  tech¬ 
niques,  including  the  group  contribution  methods  of  Lydersen  and  of  Jo- 
back.  The  solubility  of  CL20  in  supercritical  C02was  estimated  using  a 
Fortran  program  developed  during  the  course  of  this  project.  Estimations 
spanned  a  reduced  temperature  range  of  1.003  to  1.21  and  a  reduced  pres¬ 
sure  range  of  1.01  to  2.06  with  respect  to  carbon  dioxide.  The  Fortran  pro¬ 
gram  was  validated  using  available  literature  data  for  the  solubility  of 
naphthalene  and  of  biphenyl  in  supercritical  CO2.  The  applicability  of  the 
estimation  techniques  employed  for  the  critical  properties  for  CL20  was 
established  using  these  same  techniques  to  estimate  the  critical  properties 
of  comparable  compounds,  including  RDX  and  HMX.  Solubility  data  for 
RDX  in  supercritical  CO2  reported  in  the  literature  were  also  used  to  estab¬ 
lish  the  validity  of  the  estimation  approach. 
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1  Introduction 

Background 

CL20,  a  high-energy  explosive  compound,  is  a  polyazapolycyclic  caged  polyni- 
tramine  (2,4,6,8,10,12-Hexanitro-  2,4,5,8,10,12-hexaazaisowurtzitane).  Of  the 
several  crystalline  phases  of  CL20,  the  8-CL20  has  highest  density  (2.044 
g/cm3)  and  thermal  stability  (Nielson  et  al.  1998).  The  combustion  and  deto¬ 
nation  characteristics  of  the  CL20  can  be  improved  if  it  is  formed  into 
nanoparticles  of  uniform  size.  Conventional  methods  to  form  ultrafine  parti¬ 
cles  with  narrow  size  distributions  generally  follow  either  size  reduction  or 
comminution  by  direct  generation  routes  such  as  dispersion  reactions,  spray 
drying,  or  controlled  recrystallization.  Most  of  these  processes  are  energy  in¬ 
tensive,  generate  hazardous  waste  streams,  and  result  in  particles  with  wide 
size  distribution  and  shapes.  A  new  promising  process  using  environmentally 
benign  compressed  gases  as  either  solvents  or  anti-solvents  is  being  investi¬ 
gated  for  applications  in  pharmaceutical  and  other  industries.  This  process 
shows  that  ultrafine  particles  can  be  developed  using  supercritical  fluids 
(SCF)  without  generating  large  quantities  of  solvent  wastes.  This  process  may 
be  termed  as  a  green  manufacturing  process  for  producing  nanocrystalline 
particles  of  8-CL20.  Information  on  various  aspects  of  fluid  dynamics,  ther¬ 
modynamics,  and  mass  transfer  including  nucleation  and  atomization,  crystal 
growth,  and  agglomeration  of  CL20  in  supercritical  carbon  dioxide  (CO2)  are 
required  to  develop  such  a  process.  This  project  was  focused  to  develop  a 
model  for  prediction  of  solubility  of  CL20  in  supercritical  CO2. 

The  solubility  of  CL20  in  supercritical  CO2  must  be  estimated  to  determine 
the  feasibility  of  using  supercritical  processing  for  production  of  explosive 
grade  CL20.  The  solubility  of  a  solute  in  a  supercritical  fluid  can  be  readily 
evaluated  using  an  equation  of  state,  but  necessary  thermodynamic  informa¬ 
tion  required  for  the  evaluation  must  be  either  available  in  the  literature  or 
predicted  using  available  techniques.  The  approach  taken  in  this  work  was  to 
ascertain  what  data  were  available  in  the  literature  for  the  compound  of  inter¬ 
est  and  for  other  similar  energetic  compounds.  The  validity  of  available  esti¬ 
mation  techniques  for  the  required  thermodynamic  properties  was  then  as¬ 
sessed  through  comparison  of  the  predictions  with  literature  data  where 
available. 
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Objective 

The  specific  objective  of  the  research  was  to  develop  a  thermodynamic  model 
to  obtain  solubility  of  8-CL20  in  supercritical  CO2.  The  model  thus  developed 
shall  be  validated  with  available  experimental  solubility  of  similar  nitroaro- 
matic  explosives  in  supercritical  CO2. 

Approach 

To  obtain  solubility  of  the  CL20  in  supercritical  CO2  and  solutions  of  CO2  and 
co-solvents  selected,  based  upon  prior  experience  of  ARDEC  and  literature 
information,  involves  estimating  critical  properties  (from  literature  or  use 
group  contribution  methods),  and  thermophysical  and  transport  properties  of 
the  compounds  involved.  The  phase  behavior  of  the  system  was  modeled  in 
FORTRAN  computer  language  with  a  cubic  equation  of  state  such  as  Peng- 
Robinson  Equation  of  State  (Peng  and  Robinson  1976)  or  any  such  equation 
of  state.  The  solubility  data  obtained  in  the  literature  were  used  to  obtain  in¬ 
teraction  parameters  of  the  equation  of  state. 

Mode  of  Technology  Transfer 

Information  contained  in  this  report  can  be  used  as  a  basis  for  developing 
green  processes  to  produce  high  density  particulate  explosives.  The  results  of 
this  research  help  in  reducing  waste  streams  generated  during  the  manufac¬ 
turing  processes.  The  results  of  this  research  will  be  made  available  to  the 
Armament  Research  and  Development  Center  (ARDEC)  and  to  the  Pollution 
Prevention  and  Environmental  Quality  Technology  (EQT)  Programs.  Results 
of  this  research  are  also  presented  at  the  24th  Army  Science  Conference. 

This  report  will  be  accessible  through  the  World  Wide  Web  (WWW)  URL: 
http :  /  /www.  cecer.army.mil 
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2  Literature  Review 

A  comprehensive  review  of  the  literature  was  conducted  to  identify  sources  of 
thermodynamic  information  for  CL20  and  similar  energetic  compounds  in¬ 
cluding  cyclotrimethylene-trinitramine  (RDX)  and  cyclotetramethylene- 
tetranitramine  (HMX).  These  other  energetic  compounds  were  included  since 
they  have  been  used  as  explosive  agents  for  a  longer  time,  so  more  informa¬ 
tion  was  likely  to  be  available.  The  available  data  for  these  compounds  could 
then  be  used  to  validate  estimation  techniques  for  the  thermodynamic  prop¬ 
erties  of  CL20.  Two  computerized  searches  were  conducted  and  yielded  dif¬ 
ferent  items  of  interest. 

Beilstein  Database 

Available  information  on  the  compounds  of  interest  was  compiled  through 
use  of  the  Chemical  Abstracts  Service  (CAS)  registry  numbers  (CL20,  CAS# 
135285-90-4;  RDX,  CAS#  121-82-4;  HMX,  CAS#  2691-41-0).  Of  particular 
interest  was  the  determination  of  what  data  were  available  in  the  literature, 
including:  the  critical  point  (temperature,  pressure,  volume,  compressibility), 
the  melting  point,  the  acentric  factor,  the  latent  heat  of  vaporization,  vapor 
pressure,  and  liquid  and/or  solid  density.  For  CL20,  the  only  relevant  data 
contained  in  the  database  were  crystal  density  (solid  density)  and  melting 
point.  No  references  were  reported  associated  with  the  critical  point,  the 
normal  boiling  point,  or  phase  change  properties  (latent  heat,  vapor  pres¬ 
sure).  For  RDX,  relevant  references  were  found  for  the  critical  temperature, 
critical  volume,  vapor  pressure,  and  melting  point.  The  values  obtained  from 
the  references  are  compiled  in  Table  1. 

For  HMX,  relevant  references  were  also  identified  for  these  properties,  ex¬ 
cepting  vapor  pressure.  The  values  obtained  from  the  references  are  compiled 
in  Table  2. 


Table  1.  Thermodynamic  information  for  RDX  (compiled  from  Beilstein  Search). 


Property 

Value 

Reference 

Critical  Temperature 

567  °C 

Maksimov  (1992) 

Critical  Volume 

442000  cm^/gmol 

Maksimov  (1992) 

Vapor  Pressure 

0.000036-0.0004  mm  Hg@ 
110.6-138.5  °C 

Edwards  (1953) 

Melting  Point 

Range:  272  to  278  °C 

Burov  etal.  (1999) 
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Table  2.  Thernnodynannic  information  for  HMX  (compiled  from  Beilstein  Search). 


Property 

Value 

Reference 

Critical  Temperature 

654  °C 

Maksimov  (1992) 

Critical  Volume 

611000  cmygmol 

Maksimov  (1992) 

Enthalpy  of  Vaporization 

158599.97  J/mol  @  234  °C 

Behrens  (1990) 

Melting  Point 

Range:  198  to  205°C 

Suri  and  Chapman  (1988) 
Burov  et  al.  (1999) 

Chemical  Abstracts  Database 

The  Chemical  Abstracts  Database  revealed  additional  information  for  the 
compounds  of  interest.  For  RDX,  two  items  of  significance  were  identified. 
Dionne  et  al.  (1986)  measured  the  vapor  pressure  of  RDX  as  a  function  of 
temperature  over  the  temperature  range  37  to  102  °C.  Morris  (1998)  meas¬ 
ured  the  solubility  of  RDX  in  dense  CO2  over  the  temperature  range  of  303  to 
353  K.  These  data  provide  a  means  to  assess  the  validity  of  the  approach  taken 
to  estimate  vapor  pressure  and  to  predict  the  solubility  in  supercritical  CO2 
for  the  compound  of  interest.  The  majority  of  other  citations  obtained 
through  this  search,  while  of  interest,  were  not  relevant  to  the  tasks  involved 
in  this  work. 
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3  Evaluation  of  Solubility  in  a  Supercritical 
Fluid 

The  basis  for  predicting  the  solubility  of  a  solute  in  a  supercritical  fluid  sol¬ 
vent  is  the  equivalence  of  fugacities  for  the  particular  solute  in  each  phase: 

=  f-  (Eq  1) 


where  the  superscript  s  represents  the  solid  phase  and  f  the  supercritical  fluid 
phase.  If  the  solubility  of  the  supercritical  fluid  in  the  solid  phase  is  assumed 
negligible,  then  the  fugacity  of  the  solute  in  the  solid  phase,  /;%  is  equal  to  the 
fugacity  of  the  pure  solute,  f  -^ .  The  fugacity  of  the  pure  solute  in  the  solute 
phase  is  evaluated  using  (Sandler  1989): 


p 

J  1  1 


/■ 


exp 


f 

•>psa 


RT 


dP 


(Eq2) 


The  molar  volume  of  the  pure  solute,  Vj®,  is  assumed  constant.  The  fugacity  of 
the  pure  solute  in  the  saturated  state,  is  combined  with  the  pressure  in 
the  denominator  to  form  the  fugacity  coefficient  of  the  pure  solute  in  the  satu¬ 
rated  state,  Integration  of  Equation  2  with  these  substitutions  yields 
(Sandler  1989): 


/.■=P“<|.“exp|^(p-P“) 


(Eqs) 


For  many  low  volatility  compounds,  the  fugacity  coefficient  for  the  pure  solute 
in  the  saturated  state  is  very  nearly  unity  due  to  the  extremely  low  vapor  pres¬ 
sure  of  these  compounds  at  ambient  conditions. 

The  fugacity  of  the  solute  in  the  supercritical  fluid  phase  is  given  by: 

/i'=y,$iP  (Eqq) 

where  yi  is  the  mole  fraction  of  solute  in  the  supercritical  fluid  phase,  also  de¬ 
fined  as  the  solubility  of  the  solute  in  the  supercritical  fluid;  |,is  the  fugacity 
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coefficient  for  the  solute  in  the  supercritical  fluid  phase;  and  P  is  the  system 
pressure. 


Combination  of  Equations  3  and  4  yields  an  expression  for  the  solubility  of 
the  solute,  yi,  in  the  supercritical  fluid  (Sandler  1989): 


^psat  Y  ^ 

1 

-Y 

1 

_ 1 

p  YJy 

RT 

(Eqs) 


The  ideal  solubility  of  the  solute,  a  function  of  temperature  and  pressure,  is 
represented  by  the  first  bracketed  term  in  the  expression.  For  the  compounds 
of  interest,  the  ideal  solubility  will  be  extremely  low  due  to  the  low  vapor 
pressure  of  the  solute  and  the  high  pressure  required  to  achieve  supercritical 
conditions.  Non-ideal  behavior  of  the  supercritical  fluid  phase  is  represented 
by  the  second  bracketed  term  in  the  expression.  The  last  term  is  the  Poynting 
factor,  which  represents  the  effect  of  pressure  on  the  solid  phase.  The  last  two 
bracketed  terms,  when  combined,  are  known  as  the  enhancement  factor.  This 
factor  represents  the  increase  in  solubility  due  to  the  solvent’s  supercritical 
state. 


The  fugacity  of  the  solute  in  the  supercritical  fluid  can  be  evaluated  using  a 
cubic  equation  of  state  such  as  the  Redlich-Kwong-Soave  or  the  Peng  Robin¬ 
son  equations.  The  Peng-Robinson  equation  is  used  in  this  work.  The  Peng- 
Robinson  equation  for  a  mixture  is  (Modell  and  Reid  1983): 


Rh  _ a^j _ 

V-b^~V(V  +  b^)  +  b^(V-b^) 


(Eq6) 


where  V  is  the  molar  volume  of  the  mixture  while  an  and  bm  are  mixture- 
dependent  parameters.  For  the  Peng-Robinson  equation  of  state,  these  mix¬ 
ture  parameters  are  evaluated  from  mixing  rules  where  the  pure  component 
analogues,  calculated  from  the  critical  properties,  the  acentric  factor  co,  and 
the  reduced  temperature  for  each  species,  are  combined  through  appropriate 
combinatorial  relationships.  The  standard  Peng-Robinson  mixing  rules  are 
used. 


N  N 

aM=ZZy.yjay  (Eqy) 

i=l  j=l 

N 

bw  =  Zy.b, 

i=l 


(Eq8) 
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with:  a.j  =  (l -  and  a;;  =  a-.  The  binary  interaction  parameter,  5^ ,  is 

specific  for  the  i-j  binary  pair.  In  these  expressions,  the  ai  and  hi  represent  the 
pure  component  parameter  values  which  are  evaluated  using: 


a,(co,TJ  =  a(TJa(co,TJ 

0.45724R"t" 
a(Te)  = - B ^ 


(Eqp) 

(Eq  10) 


a{(i),T)  = 


1  +  (0.37464  +  i.54226(y  -  0. 269920)^  )ii-  ^ 

l  +  /c(l-7^ 


(Eqii) 


b; 


0.07780RT^ 


(Eq  12) 


The  Peng-Robinson  equation  is  often  written  in  terms  of  the  compressibility, 
Z  =  (PV/RT): 

+Z"(B-1)  +  z(-3B"  -2B  +  a)  +  (b^  +B"  -Ab)  =  0  (Eq  13) 


where  A  and  B  are  defined  as: 


A  = 


(RT)" 


(Eq  14) 


B  = 


bNiP 

RT 


(Eq  15) 


The  use  of  this  equation  of  state  for  evaluation  of  the  fiigacity  coefficient  of 
the  solute  in  a  supercritical  fluid  is  well-documented  in  the  literature 
(McHugh  and  Paulaitis  1980).  Following  the  notation  in  Modell  and  Reid 
(1983),  the  fiigacity  coefficient  for  species  i  in  solution  is  found  through: 


^(Z-i)-ln(Z-B) 


4  ( ^Hyk^ik 

^  k 

2V2B 


A 

b 


M 


In 


Z  +  b(i-V2) 

Z  +  B(l  +  ^/2) 


(Eqi6) 
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The  prediction  of  solubility  by  this  set  of  equations  requires  an  iterative 
approach.  In  general,  the  temperature  and  pressure  are  known,  as  well  as  an 
estimate  of  the  solubility,  yi.  If  experimental  data  are  not  available  for  use  as 
the  initial  estimate,  then  an  estimate  is  calculated  by  assuming  that  ideal 

behavior  is  valid  (i.e.,  y i  “  The  estimated  yi  is  then  used  with  the 

given  temperature  to  evaluate  the  pure  component  and  mixture  parameters. 
The  compressibility  is  then  determined  through  Equation  13  using  the  given  T 
and  P  with  the  mixture  parameters  through  solution  of  the  cubic  equation  of 
state.  In  this  implementation,  the  cubic  equation  is  solved  analytically.  The 
compressibility  is  substituted  into  the  expression  for  the  fugacity  coefficient. 
Equation  16,  along  with  the  necessary  pure  component  and  mixture 
parameters.  The  predicted  solubility  is  then  calculated  through  Equation  5. 
The  iterative  procedure  is  continued  until  the  predicted  solubility  at  the  end 
of  an  iteration  is  equal  to  the  estimate  at  the  start  of  the  iteration  to  within 
some  prescribed  tolerance. 
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4  Estimation  of  Criticai  Properties  and  Other 
Required  Thermodynamic  information 


Estimation  of  the  critical  properties  and  other  required  thermodynamic  in¬ 
formation  was  undertaken  for  the  energetic  compounds  RDX,  HMX,  and 
CL20.  Since  a  limited  amount  of  pertinent  data  was  available  for  RDX  and 
HMX  (i.e.,  critical  temperature,  critical  volume,  vapor  pressure),  the  estima¬ 
tion  of  properties  for  these  compounds  could  serve  as  a  baseline  to  establish 
the  validity  of  the  estimation  procedures  used  and  thus  provide  for  greater 
confidence  in  the  estimates  obtained  for  CL20.  Although  CL20  (Figure  1)  is  a 
much  more  complex  molecule  when  compared  with  RDX  (Figure  2)  or  HMX 
(Figure  3),  all  three  molecules  contain  -NO2  groups,  individually  attached  to  a 
ring  >N-  group. 
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Figure  1.  Molecular  structure  of  CL20. 
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Figure  2.  Molecular  structure  of  RDX. 
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Figure  3.  Molecular  structure  of  HMX. 


Group  contribution  methods  were  used  to  estimate  the  necessary  values.  Ta¬ 
ble  3  shows  the  group  decomposition  for  these  three  compounds. 


Table  3.  Group  decomposition  of  energetic  molecules. 


Group  ID 

CL20 

RDX 

HMX 

Number  of  Type 

Number  of  Type 

Number  of  Type 

>N-  (ring) 

6 

3 

4 

-NO2 

6 

3 

4 

CH  (ring) 

6 

CH2  (ring) 

3 

4 

The  method  of  Stein  and  Brown  (1994)  was  used  to  estimate  the  normal  boil¬ 
ing  point  for  each  compound.  This  method  is  an  extension  of  Joback’s  method 
(Joback  and  Reid  1987)  for  estimation  of  the  normal  boiling  point.  Stein  and 
Brown’s  extension  specifically  involves  inclusion  of  additional  groups  and  re¬ 
vision  of  Joback  and  Reid’s  groups  into  smaller  subdivisions.  Of  interest  in 
this  work  is  the  inclusion  of  the  >N-  (ring)  group.  For  this  group,  the  contri¬ 
butions  to  critical  temperature  and  pressure  were  identical  to  those  for  a  >N- 
(non-ring)  group,  but  were  different  for  the  normal  boiling  point.  Stein  and 
Brown  use  the  original  relation  of  Joback  and  Reid,  given  as  Equation  17: 

T^=  198.2 +  ynig.  (Eqi7) 


An  additional  refinement  to  the  estimate  is  also  given  by  Stein  and  Brown  and 
was  used  in  this  work.  For  estimated  normal  boiling  points  for  700  K  and  be¬ 
low,  the  corrected  boiling  point  is  given  by: 


Tb  (corr)  =  T,  -  94.84  +  O.SSTTT^  -  0.00077054^" 


(Eqi8) 
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Above  700  K,  the  corrected  boiling  point  is  given  by: 

T,  (corr)  =  T,  +  282.7  -  0.52091^  (Eq  19) 

Estimated  normal  boiling  points  calculated  in  this  work  for  the  compounds  of 
interest  are  shown  in  Table  4. 


Table  4.  Estimated  normal  boiling  points. 


Compound 

Tb  Estimate 
(K) 

Tb(corr)  Estimate 
(K) 

RDX 

717.8 

626.6 

HMX 

891.0 

709.6 

CL20 

1208.7 

861.8 

The  estimation  of  critical  temperature  was  accomplished  using  Joback’s 
method  (Joback  and  Reid  1987)  as  implemented  in  the  Cranium  software 
(Molecular  Knowledge  Systems,  Inc.,  Bedford,  NH).  This  estimation  tech¬ 
nique  requires  the  use  of  the  normal  boiling  point;  thus,  the  values  shown  in 
Table  4  were  used.  Table  5  shows  estimates  for  the  critical  temperature  for  the 
compounds  of  interest. 

The  critical  temperature  estimated  for  RDX,  842.7  K,  using  the  Tb(corr)  esti¬ 
mate  is  in  excellent  agreement  with  the  experimental  value  of  567  °C  (840  K) 
reported  by  Maksimov  (1992).  Similarly,  the  critical  temperature  estimated 
for  HMX,  913.5  K,  using  the  Tb(corr)  estimate  is  in  good  agreement  with  the 
experimental  value  of  654  °C  (927  K)  reported  by  Maksimov  (1992). 


Table  5.  Estimated  critical  temperatures. 


Compound 

Critical  Temperature,  Tc 
Estimated  using  Tb 
(K) 

Critical  Temperature,  To 
Estimated  using  Tb(corr) 

(K) 

RDX 

965.4 

842.7 

HMX 

1147.1 

913.5 

CL20 

1483.6 

1057.8 
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The  estimation  of  critical  pressure  was  also  accomplished  using  Joback’s 
method  (Joback  and  Reid  1987)  as  implemented  in  the  Cranium  program.  No 
literature  data  for  the  critical  pressure  were  identified  for  the  compounds  of 
interest.  Table  6  shows  estimates  for  the  critical  temperature  for  the  com¬ 
pounds  of  interest. 


Table  6.  Estimated  critical  pressure. 


Compound 

Critical  Pressure,  Pc 
(bar) 

RDX 

58.0 

HMX 

53.0 

CL20 

48.9 

In  addition  to  the  critical  properties  (temperature  and  pressure),  the  acentric 
factor,  defined  as  vapor  pressure  of  the  compound  at  a  reduced  pressure  of 
0.7,  is  used  in  describing  the  temperature  dependence  of  the  attraction  term 
in  the  Peng-Robinson  equation  of  state.  Thus,  estimation  of  the  acentric  fac¬ 
tor  was  also  required  for  the  compounds  of  interest.  Poling  et  al.  (2001)  rec¬ 
ommend  the  estimation  of  the  acentric  factor  using  the  three  parameter  Pitzer 
expansion: 

ln(P,p  /PJ  =  +cof<'>  +co"f<">  (Eq  20) 


Neglecting  the  second  order  term,  this  can  be  rearranged,  as  shown  in  Equa¬ 
tion  21  to  provide  an  estimate  of  the  acentric  factor  using  the  critical  point 
and  the  normal  boiling  point. 


CO  =  - 


ln(P^/1.01325)  +  f<°^(Tb,) 
f^CTJ 


(Eq  21) 


where  Tbr  is  the  reduced  normal  boiling  point  and  pressure  is  in  bar.  The 
functions,  Eo)(Tbr),  ^^^(Tbr),  and  E2)(Tbr)  are  the  temperature  dependent 
correlations  developed  by  Ambrose  and  Walton  (1989)  for  the  three 
parameter  corresponding  states  method  for  prediction  of  vapor  pressure. 
These  expressions  are  shown  in  Equations  22  to  24  (Poling  et  al.  2000). 

^(0) ^  _  -5.97616t  +  1.29874t^-^  -0.60394t"-^  -1 .068411^ 


(Eq  22) 
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-5.03365t  +  1.11505t^-^  -5.41217t^^  -7.46628t^ 

T. 

-0.64771t  +  2.41539t‘^  -4.26979t"-^  +3.25259t^ 

T 

r 


(Eq  23) 

(Eq  24) 


In  these  expressions,  t  is  defined  as  (i-Tr).  For  larger  molecules,  the  inclusion 
of  the  second  order  term  can  be  important.  Thus,  inclusion  of  the  second  or¬ 
der  term  gives  rise  to  a  quadratic  equation  that  must  be  solved.  This  quadratic 
equation  is  in  the  form: 

-ln(l. 01325 /PJ  =  0  (Eq  25) 


The  acentric  factor  was  estimated  using  the  values  of  Tb(corr)  shown  in  Table 
4,  and  values  of  Tc  from  Table  5.  Estimated  values  for  each  of  the  compounds 
of  interest  are  shown  in  Table  7.  These  values  are  compared  to  a  two-point  ex¬ 
trapolation  of  the  vapor  pressure  curve  to  obtain  an  estimate  of  the  acentric 
factor.  Assuming  a  linear  relationship  between  logio(Pr)  and  i/Tr  gives  rise  to 
the  expression  shown  in  Equation  26. 

(Eq  26) 

where  B  is  given  by: 


logio(Pr)  =  -^ 


B  = 


10glo(Pc)-loglo(Pnbp) 


1 


^  nbp 


1 

T. 


(Eq  27) 


Using  the  estimated  normal  boiling  point  and  the  critical  temperature  and 
pressure,  the  reduced  vapor  pressure  at  a  reduced  temperature  of  Tr  =  0.7  was 
evaluated  and  then  used  in  the  definition  of  the  acentric  factor.  Estimated 
values  for  the  acentric  factor  obtained  using  this  method  are  also  provided  in 
Table  7. 
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Table  7.  Estimated  acentric  factor. 


Compound 

Acentric  Factor 
(Estimated  Using  Cranium) 

Acentric  Factor 
(Estimated  Using 

Eq.  22-25) 

Acentric  Factor 
(Estimated  Using 

Eq.  26  and  27) 

RDX 

1.25 

1.22 

1.19 

HMX 

1.57 

1.66 

1.80 

CL20 

2.40 

2.38 

2.18 

The  estimation  of  vapor  pressure  was  also  undertaken  for  the  model  com¬ 
pounds.  Vapor  pressure  data  were  located  for  RDX,  but  data  for  HMX  and 
CL20  were  not  identified  through  the  literature  search.  Dionne  et  al.  (1986) 
report  the  vapor  pressure  of  RDX  as 

log,„(P.,)  =  22.50--^  (Eq28) 

where  Pvp  is  in  parts  per  trillion  (v/v)  and  T  is  in  K.  This  correlation  was  based 
on  their  measurements  over  the  temperature  range  of  37  to  102  °C  combined 
with  literature  data  from  three  other  sources.  A  second  expression  for  the  va¬ 
por  pressure  of  RDX  was  (U.S.  Dept,  of  Labor  2003): 

log,o(Pvp)  =  14.14-^  (Eq29) 

where  Pvp  is  in  Pa  and  T  is  in  K.  These  two  expressions  provide  roughly  the 
same  estimates  for  the  vapor  pressure  of  RDX  over  the  temperature  range  of 
interest. 

Poling  et  al.  (2001)  recommend  the  use  of  the  corresponding  states  method 
with  the  Ambrose-Walton  expressions  for  the  estimation  of  vapor  pressure. 
The  relevant  expressions  are: 

ln(Pvpj.)  =  +cof*'^ +cof^^^  (Eq30) 

where  Pvpr  is  the  reduced  vapor  pressure,  is  given  by  Equation  22,  by 
Equation  23,  and  by  Equation  24.  Vapor  pressures  were  estimated  accord¬ 
ing  to  Equation  30  using  the  values  for  the  acentric  factor  given  in  Table  7. 
Vapor  pressures  for  RDX  estimated  using  this  technique  differed  from  the 
available  literature  data  by  a  factor  of  io3.  Vapor  pressures  were  also  esti¬ 
mated  according  to  Equation  26.  The  estimated  vapor  pressure  is  compared 
to  the  available  literature  data  for  RDX  in  Eigure  4.  The  estimates  are  in  good 
agreement  with  the  literature  data  over  the  range  of  25  °C  to  approximately 
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6o  °C.  Above  6o  °C,  the  estimated  vapor  pressure  and  the  literature  data  di¬ 
verge,  with  the  estimated  vapor  pressure  being  almost  half  that  reported  in 
the  literature  at  8o  °C.  For  prediction  of  solubility  in  supercritical  CO2,  it  is 
anticipated  that  the  temperature  range  of  interest  is  between  approximately 
35  and  60  °C.  Thus,  Equation  26  was  used  to  estimate  the  vapor  pressure  of 
CL20.  These  estimated  vapor  pressures  were  then  used  in  the  evaluation  of 
the  solubility  of  CL20  in  supercritical  CO2  and  are  plotted  in  Figure  5. 
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Figure  4.  Comparison  of  estimated  RDX  vapor  pressure  with  iiterature  data. 
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Figure  5.  Estimated  vapor  pressure  for  CL20. 
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5  Validation  of  the  Deveioped  Code 

The  Fortran  code  was  validated  using  the  solubility  of  naphthalene  and  of  bi¬ 
phenyl  in  supercritical  CO2.  Data  for  these  systems  were  measured  by 
McHugh  and  Paulaitis  (1980).  For  naphthalene  in  supercritical  CO2,  the  value 
of  the  interaction  parameter,  812,  was  set  to  0.103,  while  for  biphenyl  in  super¬ 
critical  CO2,  a  value  of  0.12  was  used.  These  values  for  the  interaction  parame¬ 
ter  were  obtained  by  Schmitt  (1984)  by  regression  of  the  data  of  Tsehanskaya 
et  al.  (1964)  for  naphthalene,  and  the  data  of  McHugh  and  Paulaitis  for  bi¬ 
phenyl.  Figure  6  provides  a  comparison  of  the  literature  data  to  the  predicted 
values  for  naphthalene,  while  Figure  7  provides  the  comparison  of  literature 
data  to  predicted  values  for  biphenyl.  The  predicted  values  are  in  fair  agree¬ 
ment  with  the  literature  data  for  naphthalene  in  the  range  of  pressures  from 
approximately  the  critical  pressure  of  CO2  (~  74  bar)  to  approximately  200 
bar.  The  predicted  solubilities  at  55  °C,  corresponding  to  the  highest  pres¬ 
sures  examined  by  McHugh  and  Paulaitis,  diverge  from  the  experimental 
data.  A  reason  for  the  discrepancy  is  that  Schmitt  (1984)  reported  different 
values  of  the  interaction  parameter  for  each  set  of  isothermal  data.  Thus,  the 
same  interaction  parameter  did  not  provide  the  best  representation  of  data 
for  differing  isotherms. 
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Figure  6.  Solubility  of  naphthalene  in  supercritical  CO2;  comparison  of  literature  data  and 

calculated  values. 

For  biphenyl  in  supercritical  CO2,  the  agreement  between  predicted  and 
measured  values  is  much  better,  which  is  most  likely  due  to  the  use  of  an 
interaction  parameter  for  the  data  set  that  had  been  obtained  through 
regression  of  the  same  data  set.  The  primary  motivation  of  using  these  well- 
known  data  sets  to  provide  validation  and  verification  of  the  developed  code 
was  that  the  performance  of  the  code  could  be  more  easily  established  since 
all  required  information  (i.e.,  critical  temperature,  critical  pressure,  acentric 
factor,  vapor  pressure,  saturated  molar  volume)  was  available  for  the  solutes 
(naphthalene  and  biphenyl)  as  well  as  the  solvent  (CO2). 


Mole  Fraction  Biphenyl  in  SC  Carbon  Dioxide,  y1 


ERDC/CERLTR-06-13 


0.018 

0.016 

0.014 

0.012 

0.01 

0.008 

0.006 

0.004 

0.002 

0 


A 


A 


A 


■  Literature,  McHugh  and  Paulaitis 
A  Calculated,  This  Work 


50 


100 


150 


200  250  300 

Pressure,  P,  (bar) 


350 


400 


450 


500 


Figure  7.  Solubility  of  biphenyl  in  supercritical  CO2  comparison  of  literature  data  and  calculated 

values  at  35.8  °C. 
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6  Solubility  of  Explosives  in  Supercritical  C02 

Comparison  of  SCF  Predictions  for  RDX  in  CO2  With  Avaiiabie  Literature 
Data 

The  solubility  of  RDX  in  subcritical  and  supercritical  carbon  dioxide  was  re¬ 
ported  by  Morris  (1998).  The  reported  values  (in  mg  RDX/g  CO2)  were  con¬ 
verted  to  mole  fractions  for  comparison  with  code  predictions.  The  critical 
temperature,  critical  pressure,  and  acentric  factor  estimated  for  this  work 
were  used  along  with  the  vapor  pressure  reported  by  Dionne  et  al.  (1986)  to 
predict  the  solubility  of  RDX.  The  initial  estimate  used  in  the  iterative  code 
for  the  solubility  was  the  experimental  value  reported  by  Morris  at  the  condi¬ 
tions  of  interest.  There  is  little  agreement  between  the  literature  data  and  the 
predictions,  as  shown  in  Figure  8.  The  literature  solubilities  are  quite  small, 
with  magnitudes  between  lo  ^  and  lo  ^,  while  the  predicted  solubilities  were 
larger  with  magnitudes  between  lo  ^  and  io-3.  This  is  most  likely  due  to  the 
vapor  pressure  expression  employed.  Error  in  the  vapor  pressure  translates 
directly  into  error  in  the  predicted  solubility,  as  evidenced  through  Equation 
5.  If  the  vapor  pressure  is  overestimated  by  a  factor  of  2,  then  the  predicted 
solubility  will  be  overestimated  to  the  same  degree.  Accurate  measurement  of 
the  vapor  pressure  of  the  explosive  compound  is  an  essential  requirement  for 
reliable  predictions  of  the  solubility  of  an  explosive  in  supercritical  CO2. 
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Figure  8.  Solubility  of  RDX  in  subcritical  and  supercritical  CO2;  comparison  of  predictions  with 

literature  data. 

SCF  Predictions  for  CL20  in  Supercriticai  CO2 

The  solubility  of  CL20  in  supercritical  CO2  was  predicted  using  the  developed 
Fortran  code.  The  estimated  critical  temperature,  critical  pressure,  and  acen¬ 
tric  factor  were  used  in  the  code  along  with  Equation  26  for  estimation  of  the 
vapor  pressure.  A  two-point  fit  (critical  point,  normal  boiling  point)  allowed 
estimation  of  the  coefficient  B  in  Equation  26.  Since  no  experimental 
measurements  were  available,  the  initial  estimate  of  the  solubility  was  set  to 
i.oe-16.  Solubility  was  predicted  over  the  temperature  range  of  305.15  to 
368.15  K  and  over  the  pressure  range  of  74  to  150  atmospheres  (atm).  The 
temperature  range  corresponds  to  a  range  of  reduced  temperature  of  1.003  to 
1.21  with  respect  to  CO2,  while  the  pressure  range  corresponds  to  a  reduced 
pressure  range  of  1.01  to  2.06.  The  estimated  solubilities  are  shown  in  Eigure 
9  (low  temperature)  and  Eigure  10  (high  temperature).  In  general,  as  the 
temperature  increases,  the  solubility  decreases,  while,  as  the  pressure 
increases,  the  solubility  increases.  The  predicted  behavior  is  markedly  dif¬ 
ferent  near  the  critical  point  of  CO2  when  compared  with  the  behavior  of 
naphthalene  in  supercritical  CO2.  An  examination  of  Eigure  6  indicates  that, 
as  the  pressure  is  increased  from  near  the  critical  point  to  higher  pressures, 
the  solubility  of  naphthalene  shows  an  inflection  point,  followed  by  a 
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significant  increase  in  the  solubility  over  a  fairly  narrow  pressure  range.  This 
increase  is  influenced  to  a  significant  degree  by  the  vapor  pressure  of  the 
solute.  For  CL20,  the  estimated  vapor  pressures  are  extremely  small,  on  the 
order  of  10-18  at  ambient  temperature,  increasing  to  10-13  at  368.15  K.  Thus, 
the  predicted  solubilities  are  also  small  (range  of  10-13  to  lo-^),  with  the 
highest  solubility  predicted  for  308.15  K  (35°C)  and  150  atm.  In  the  region  of 
temperature  nearest  the  critical  point  of  CO2,  the  influence  of  the  supercritical 
fluid  is  stronger.  It  would  be  most  desirable  to  operate  in  this  region  of  35- 
50  °C  if  one  is  interested  in  maximizing  the  solubility  of  the  solute,  CL20,  in 
the  CO2. 


Pressure,  P,  (bar) 


Figure  9.  Predicted  solubility  of  CL20  in  supercritical  CO2  -  low  temperature  range. 
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7  Conclusion 

A  comprehensive  search  of  the  literature  revealed  limited  data  that  could  be 
used  to  validate  the  approach  taken  in  this  work  for  estimation  of  necessary 
thermodynamic  properties.  The  estimation  of  normal  boiling  point,  critical 
temperature,  critical  pressure,  acentric  factor,  and  vapor  pressure  for  CL20 
was  accomplished  using  available  group  contribution  techniques  and  property 
definitions.  While  CL20  is  a  very  complex  species,  comparison  of  estimated 
critical  temperature  for  RDX  and  HMX  to  values  reported  in  the  literature 
provides  a  measure  of  confidence  in  the  estimated  values  for  CL20.  All  three 
species  are  composed  of  similar  molecular  groups,  the  primary  difference  be¬ 
ing  the  presence  of  -CH  (ring)  groups  in  CL20  and  -CH2  (ring)  groups  in  RDX 
and  HMX. 

A  Fortran  code  was  developed  for  the  estimation  of  the  solubility  of  a  solute  in 
supercritical  CO2.  Performance  of  the  code  was  validated  using  the  data  of 
McHugh  and  Paulaitis  (1980)  for  naphthalene  and  biphenyl  solubilities  in  su¬ 
percritical  CO2.  The  prediction  of  solubility  for  the  explosive  RDX  in  super¬ 
critical  C02was  performed  and  the  predicted  values  compared  with  available 
literature  data.  The  lack  of  agreement  between  the  predicted  and  literature 
data  is  most  likely  due  to  the  uncertainty  associated  with  the  vapor  pressure 
of  RDX. 

Solubilities  for  CL20  in  supercritical  CO2,  over  a  reduced  temperature  range 
of  1.003  to  1.21  and  a  reduced  pressure  range  of  1.01  to  2.06,  were  predicted 
using  the  developed  code.  These  predictions  lead  to  two  primary  conclusions. 
First,  the  solubility  estimates  are  very  strongly  influenced  by  the  vapor  pres¬ 
sure.  Experimental  measurement  of  the  vapor  pressure  of  CL20  should  be 
performed.  Second,  the  conditions  that  appear  to  be  most  viable  for  maximiz¬ 
ing  the  solubility  of  CL20  in  supercritical  CO2  are  the  lowest  supercritical  tem¬ 
perature  examined  (35  °C)  and  the  highest  pressure  (150  atm). 
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Appendix  A;  Validation  of  FORTRAN  Code 


Table  Al.  Calculated  solubility  of  naphthalene  In  supercritical  CO2. 


Pressure 

(atm) 

Solubility 

T  =  308.15  K 

Pressure 

(atm) 

Solubility 

T  =  328.15  K 

85.7 

0.0047556 

81.1 

0.000666 

96.9 

0.0071149 

91.1 

0.001136 

105.1 

0.0084182 

101 

0.002154 

131.3 

0.0116050 

108.3 

0.003600 

166.8 

0.0147523 

120.1 

0.007367 

196.9 

0.0168408 

131.5 

0.011752 

219.5 

0.0181707 

141.8 

0.015812 

239.3 

0.0192039 

158.5 

0.022442 

252 

0.0198117 

169.4 

0.026870 

172.6 

0.028199 

187.1 

0.034487 

207.6 

0.044719 

220.6 

0.053024 

234.7 

0.068143 

Table  A2.  Calculated  solubility  of  biphenyl  in  supercritical  CO2. 


Pressure  (atm) 

Solubility 

T  =  308.95  K 

104.6 

0.007263 

126.2 

0.009996 

165.1 

0.013270 

201.7 

0.015190 

239.1 

0.016361 

297.3 

0.017091 

353.4 

0.016984 

435.4 

0.016064 
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Appendix  B;  Calculated  Solubilities  for  Explosive 
Compounds  in  Supercritical  CO2 


Table  Bl.  Calculated  solubility  of  RDX  In  supercritical  CO2. 


Pressure 

(atm) 

Pressure 

(bar) 

305.15  K 

308.15  K 

318.15  K 

328.15  K 

338.15  K 

348.15  K 

358.15  K 

368.15  K 

74 

74.9805 

7.84e-06 

1.80e-06 

9.19e-07 

9.89e-07 

1.30e-06 

1.88e-05 

2.82e-06 

4.32e-06 

75 

75.9938 

6.03e-05 

2.42e-06 

1.03e-06 

1.07e-06 

1.38e-06 

1.97e-06 

2.93e-06 

4.47e-06 

80 

81.0600 

3.21e-04 

3.91e-05 

1.98e-06 

1.63e-06 

1.89e-06 

2.52e-06 

3.59e-06 

5.28e-06 

90 

91.1925 

1.02e-03 

4.35e-04 

1.13e-05 

4.27e-06 

3.76e-06 

4.29e-06 

5.52e-06 

7.54e-06 

100 

101.3250 

2.15e-03 

1.08e-03 

8.18e-05 

1.33e-05 

8.08e-06 

7.64e-06 

8.74e-06 

l.lOe-05 

110 

111.4575 

4.41e-03 

2.08e-03 

2.82e-04 

4.41e-05 

1.84e-05 

1.41e-05 

1.42e-05 

1.64e-05 

120 

121.5900 

3.83e-03 

6.18e-04 

1.24e-04 

4.22e-05 

2.65e-05 

2.33e-05 

2.46e-05 

130 

131.7225 

l.lle-03 

2.73e-04 

9.09e-05 

4.94e-05 

3.85e-05 

3.71e-05 

140 

141.8550 

1.80e-03 

5.03e-04 

1.77e-04 

8.91e-05 

6.29e-05 

5.57e-05 

150 

151.9875 

2.78e-03 

8.25e-04 

3.10e-04 

1.52e-04 

l.OOe-04 

8.30e-05 

Table  B2.  Calculated  solubility  of  CL20  in  supercritical  CO2. 


Pressure 

(atm) 

Pressure 

(bar) 

305.15  K 

308.15  K 

318.15  K 

328.15  K 

338.15  K 

348.15  K 

358.15  K 

368.15  K 

74 

74.9805 

6.27e-12 

4.07e-13 

1.16e-13 

1.33e-13 

2.22e-13 

4.36e-13 

9.24e-13 

2.03e-12 

75 

75.9938 

2.43e-10 

7.15e-13 

1.46e-13 

1.56e-13 

2.51e-13 

4.81e-13 

l.OOe-12 

2.18e-12 

80 

81.0600 

5.55e-09 

1.27e-10 

5.09e-13 

3.53e-13 

4.67e-13 

7.94e-13 

1.53e-12 

3.12e-12 

90 

91.1925 

4.11e-08 

1.04e-08 

1.41e-ll 

2.29e-12 

1.80e-12 

2.31e-12 

3.68e-12 

6.56e-12 

100 

101.3250 

1.28e-07 

4.92e-08 

5.83e-10 

2.03e-ll 

8.00e-12 

7.23e-12 

9.28e-12 

1.43e-ll 

110 

111.4575 

2.93e-07 

1.36e-07 

5.76e-09 

1.99e-10 

3.93e-ll 

2.40e-ll 

2.43e-ll 

3.17e-ll 

120 

121.5900 

5.66e-07 

2.94e-07 

2.38e-08 

1.40e-09 

1.94e-10 

8.21e-ll 

6.50e-ll 

7.17e-ll 

130 

131.7225 

9.83e-07 

5.51e-07 

6.54e-08 

6.17e-09 

8.45e-10 

2.75e-10 

1.74e-10 

1.63e-10 

140 

141.8550 

1.58e-06 

9.38e-07 

1.44e-07 

1.91e-08 

3.00e-09 

8.59e-10 

4.54e-10 

3.65e-10 

150 

151.9875 

2.41e-06 

1.49e-06 

2.76e-07 

4.65e-08 

8.66e-09 

2.42e-09 

1.13e-09 

8.00e-10 
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Appendix  C;  Developed  FORTRAN  Code  for 
Estimation  of  Solubility  of  CL20  in 
Supercritical  CO2 

Instructions  for  Use  of  Executable  File 

Input  data  are  provided  in  the  SCFIN.TXT  file,  which  should  exist  in  the  same 
directory  as  the  executable  code.  The  format  for  the  input  data  file  is: 

Card  1:  Columns  1-3, 13  field,  Number  of  data  points  (NDAT)  to  be  evaluated 
(maximum  value  for  NDAT  is  100) 

Card  2:  NDAT+i:  Columns  1-10,  Temperature  in  Kelvin  (F10.5  field) 

Columns  11-20,  Pressure  in  Atmospheres  (F10.5  field) 


Sample  Input  File 


10 


368 . 15 

74 . 0 

368 . 15 

75.0 

368 . 15 

80.0 

368 . 15 

90.0 

368 . 15 

100.0 

368 . 15 

110.0 

368 . 15 

120.0 

368 . 15 

130.0 

368 . 15 

140.0 

368 . 15 

150.0 

Output  from  the  program  is  in  tabular  form.  Temperature  is  in  K,  pressure  is 
in  bar,  estimated  solubility  of  CL20  is  given  as  mole  fraction  CL20  in  the 
CL20-CO2  mixture. 
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Code  for  Executable  File 

c 

C  PROGRAM  TO  EVALUATE  THE  SOLUBILITY  OF  A  SOLUTE  IN  A  SUPERCRITICAL 

C  FLUID. 

C 

IMPLICIT  REAL*8 (A-H, 0-Z) 

COMMON  /PARAM/  AIC , BIC , A2C , B2C 

COMMON  /PARAMI/  TCI , PCI , OMEGAl , TC2 , PC2 , 0MEGA2 , DEL12 
COMMON  /PARAM2/  A1 , A2 , A12 , AMIX, B1 , B2 , BMIX, CAPAMIX, CAPBMIX 
COMMON  /PARAM3/  CAPAl , CAPA2 , CAPBl , CAPB2 
REAL* 8  TEXP (100) , PEXP (100) , Y2FIN (100) 

CHARACTER*10  INFILE, OUTFIL 
C 

INFILE=" SCFIN.TXT" 

OUTFIL="SCFOUT . TXT" 

OPEN (UNIT=5, FILE=INFILE, ACTION="READ") 

OPEN (UNIT=6, FILE=OUTFIL, ACTION="WRITE") 

C 

C  EXPERIMENTAL  DATA  ARE  READ  FROM  FILE  SCFIN.TXT 

C 

READ (5, 8000)  NDAT 
8000  FORMAT (13) 

WRITE (6, 9010) 

9010  FORMAT ( '  ',2X, 'INPUT  DATA') 


DO  10  J=1,NDAT 
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READ (5, 8001)  TEXP ( J) , PEXP ( J) 

8001  FORMAT (2F10 . 5) 

PEXP ( J) =PEXP ( J) *1 . 01325D0 
WRITE (6, 8002)  J, TEXP (J) , PEXP (J) 

8002  FORMAT ( '  ' , 2X, 13, 2X, ' T  = ' , D15 . 6, ' K ' , 2X, ' P  = ' , D15 . 6, ' BAR' ) 

10  CONTINUE 

C 

C  PROPERTIES  OF  THE  SUPERCRITICAL  FLUID  ARE  NEEDED  TO  EVALUATE  EOS 

C  PARAMETERS  -  COMPONENT  1  IS  THE  SOLVENT;  COMPONENT  2  IS  THE  SOLUTE 

C 

TC1=304 . 12D0 

PC1=73 . 74D0 

OMEGA1=0.225DO 

TC2=1058.0d0 

PC2=48 . 9d0 

OMEGA2=2 . 18d0 

DEL12=0 . ODO 

C 

C  UNIVERSAL  GAS  CONSTANT  DEFINED 

C 

RG=83 . 14D0 

C 

C  EOS  PARAMETERS  MUST  BE  EVALUATED 


C 
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A1C=0 . 4572 4D0*RG*RG*TC1*TC1/ PCI 
B1C=0 . 07780D0*RG*TC1/PC1 
A2C=0 . 45724D0*RG*RG*TC2*TC2/PC2 
B2C=0 . 07780D0*RG*TC2/PC2 

C 

C  WRITE(6,200)  A1C,B1C 

C  200  FORMATC  ',2X,  'AlC  =  '  ,  F15 . 6,  2X,  '  BIG  =',F15.6) 

C 

C  WRITE(6,201)  A2C,B2C 

C  201  FORMATC  ',2X,  'A2C  =  '  ,  F15 . 6,  2X,  '  B2C  =',F15.6) 

C 

C  MAIN  LOOP  IN  THE  PROGRAM  -  THIS  LOOP  IS  COMPLETED  FOR  EACH  OF  THE 

C  NDAT  POINTS. 

C 

DO  1000  J=1,NDAT 
Y2=l .D-16 

Yl=l .D0-Y2 
T=TEXP ( J) 

P=PEXP ( J) 

C 

C  SUBROUTINE  EOSROOT  EVALUATES  THE  COMPRESSIBILITY  OF  THE  FLUID 

C  USING  THE  PENG-ROBINSON  EQUATION  OF  STATE,  GIVEN  A  VALUE  FOR 

C  TEMPERATURE  AND  PRESSURE 

C 
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ITER=1 

C  WRITE (6, 8007)  J 

C8007  FORMAT ( '  ','J  =',I3) 

20  CALL  EOSROOT (T, P, Yl, Y2, Z, DISC) 

IF(DISC.LE. (O.DO) )  THEN 
GO  TO  1000 

ENDIF 

C  WRITE (6, 8005) 

C8005  FORMAT ( '  2X, ' FINISHED  EOSROOT') 

CALL  FUGEVAL (T, P, Yl, Y2, Z, PHI2H) 

C  WRITE (6, 8006) 

C8006  FORMAT ( '  2X, ' FINISHED  FUGEVAL') 

CALL  SOLUEST ( T , P, Yl , Y2 , Z , PHI2H, Y2EST) 

DIFF=DABS ( (Y2-Y2EST) /Y2EST) 

IF (DIFF.LE. (1 .D-16) )  THEN 
C  WRITE (6, 8003)  T,P,Y1,Y2 

C8003  FORMATC  ',2X,  'T  =  '  ,  D15 . 6,  2X,  '  P  =  '  ,  D15 . 6,  2X,  '  Yl  =  '  ,  D15 . 6,  2X,  '  Y2=  ' 
C  $,D15.6) 

C  WRITE (6, 8004)  Z , PHI2H, Y2EST 

C8004  FORMATC  ',2X,  'Z  =  '  ,  D15 . 6,  2X,  '  PHI2H  =  '  ,  D15 . 6, 2X,  '  Y2EST  =',D15.6) 
Y2FIN ( J) =Y2EST 
GO  TO  1000 

ELSE 


ITER=ITER+1 
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C  WRITE (6, 8003)  T,P,Y1,Y2 

C  WRITE (6, 8004)  Z , PHI2H, Y2EST 

Y2=Y2EST 

Yl=l-Y2 

GO  TO  20 

ENDIF 

1000  CONTINUE 

C 

C  OUTPUT  OF  FINAL  RESULTS 

C 


WRITE (6, 2003) 

2003  FORMAT ( '  ') 

WRITE (6, 2003) 

WRITE (6, 2003) 

WRITE (6, 2001) 

2001  FORMAT ( '  ' , 3X, '  J  ' , 4X, ' Temperature  (K) ' , 2X, ' Pressure  (bar) ' , 2X, 

$' Solubility  (mol  CL20/mol  mixture) ') 


WRITE (6, 2003) 

DO  2000  J=1,NDAT 

WRITE (6, 2002)  J,  TEXP ( J) , PEXP ( J) , Y2FIN ( J) 
2002  FORMAT ( '  ' , 4X, 12 , 4X, F15 . 2 , 2X, F15 . 2 , 2X, D2 1 . 6) 

2000  CONTINUE 

STOP 


END 
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CAPA2 , CAPBl , CAPB2 
C 

C  IN  THE  SUPERCRITICAL  REGION,  ONLY  A  SINGLE  ROOT  SHOULD  BE  PRESENT 

C  FOR  THE  CUBIC  EQUATION  OF  STATE.  THE  ROOT  IS  DETERMINED  USING  THE 

C  THE  STANDARD  ANALYTIC  FORM  OF  A  THIRD  ORDER  POLYNOMIAL 

C 

C  PARAMETER  A1  =  AIC*ALPHA (OMEGAl , TRl )  IS  EVALUATED 

C 

RG=83 . I4D0 
TR1=T/TC1 

EKAPPAI=0 . 37464D0+I . 54226DO*OMEGA1-0 . 26992DO*OMEGAI*OMEGA1 
TERM1=I .DO-DSQRT (TRl) 

ALPHA1= (1 .D0+EKAPPAI*TERM1) **2 
A1=AIC*ALPHA1 

B1=BIC 

C 

C  PARAMETER  A2  =  A2C*ALPHA (OMEGA2 , TR2 )  IS  EVALUATED) 

C 

TR2=T/TC2 

EKAPPA2=0 . 37464D0+I . 54226D0*OMEGA2-0 . 26992DO*OMEGA2*OMEGA2 
TERM2=I .DO-DSQRT (TR2) 

ALPHA2= (1 .D0+EKAPPA2*TERM2) **2 
A2=A2C*ALPHA2 


B2=B2C 
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C 

C  WRITE (6, 300)  ALPHAl , ALPHA2 

C  300  FORMAT ( '  'ALPHAl  = ' , FI 5 . 6 , 2X, ' ALPHA2  =',F15.6) 

C  WRITE(6,301)  A1,B1 

C  301  FORMAT ( '  ' , 2X, D15 . 6, 2X, D15 . 6) 

C  WRITE(6,301)  A2,B2 

C 

CAPA1=A1*P/ ( (RG*T) **2) 

CAPB1= (B1*P) / (RG*T) 

CAPA2=A2*P/ ( (RG*T) **2) 

CAPB2= (B2*P) / (RG*T) 

C  WRITE (6, 302)  CAPA1,CAPB1 

C  302  FORMAT ( '  ' , 2X, D15 . 6, 2X, D15 . 6) 

C  WRITE (6, 302)  CAPA2,CAPB2 

C 

C  MIXTURE  PARAMETERS  ARE  EVALUATED  USING  MIXING  RULES 

C 

A12= (1 .D0-DEL12) * (DSQRT (A1*A2) ) 

AMIX=Y1*Y1*A1+2*Y1*Y2*A12+Y2*Y2*A2 

BMIX=Y1*B1+Y2*B2 
CAPAMIX=AMIX*P/ ( (RG*T) **2) 

CAPBMIX= (BMIX*P) / (RG*T) 

C  WRITE (6, 303)  A12 

C  303  FORMATC  ',2X,  'A12  =',D15.6) 
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C  WRITE (6, 304)  AMIX,BMIX 

C  304  FORMAT ( '  ',2X, 'AMIX  = ' , D15 . 6, 2X, ' BMIX  =',D15.6) 

C  WRITE (6, 305)  CAPAMIX, CAPBMIX 

C  305  FORMAT ( '  2X, ' CAPAMIX  =', DI5 . 6, 2X, ' CAPBMIX  =',DI5.6) 

C 

C  CUBIC  FORM  OF  PENG-ROBINSON  EQUATION  OF  STATE  IS  SOLVED  FOR  THE 

C  COMPRESSIBILITY  USING  ANALYTIC  SOLUTION  OF  THE  CUBIC  EQUATION 

C 

10  COEFl=CAPBMIX-l .DO 

COEF2=-3 . D0*CAPBMIX*CAPBMIX-2 . DO *CAPBMIX+CAPAMIX 

COEF3=CAPBMIX*CAPBMIX*CAPBMIX+CAPBMIX*CAPBMIX-CAPAMIX*CAPBMIX 
QP= (3 .DO*COEF2-COEF1*COEF1) /9 .DO 

RP= (9.DO*COEFl*COEF2-27 .DO*COEF3-2*COEF1*COEF1*COEF1) /54 .DO 
DISC= (QP**3) + (RP**2) 

C  WRITE(6,309)  COEFl, COEF2, COEF3, QP, RP, DISC 

C  309  FORMATC  '  ,  6  (2X,  FIO  .  5)  ) 

C  IF(DISC.LE.O)  THEN 

C  WRITE(6,308) 

C  308  FORMATC  2X, ' DISCRIMINANT  IS  LESS  THAN  OR  EQUAL  TQ  ZERQ ' ) 

C  GQ  TO  12 

C  ENDIF 

DISCRT=DSQRT (DISC) 

ARG1=RP+DISCRT 


ARG2=RP-DISCRT 
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IF(ARG1.LT. (O.DO) )  THEN 
ARG1=DABS (ARGl) 

SP=-1.D0* (ARGl** (1.D0/3.D0) ) 

ELSE 

SP=(ARG1) ** (1.D0/3.D0) 

ENDIF 

IF (ARG2.lt. (O.DO) )  THEN 
ARG2=DABS (ARG2) 

TP=-1.D0* (ARG2** (1.D0/3.D0) ) 

ELSE 

TP=(ARG2) ** (1.D0/3.D0) 

ENDIF 

Z=SP+TP- (COEF1/3.DO) 

C  WRITE(6,310)  DISCRT, SP, TP, Z 

C  310  FORMATC  '  ,  4  (2X,  D15 . 7  )  ) 

12  RETURN 

END 

C 

-k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k 

C 

SUBROUTINE  FUGEVAL (T, P, Yl, Y2, Z, PHI2H) 

IMPLICIT  REAL*8 (A-H, 0-Z) 

COMMON  /PARAM/  AlC , BIG , A2C , B2C 


COMMON  /PARAMl/  TCI , PCI , OMEGAl , TC2 , PC2 , OMEGA2 , DEL12 
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COMMON  /PARAM2/  A1 , A2 , A12 , AMIX, B1 , B2 , BMIX, CAPAMIX, CAPBMIX 
COMMON  /PARAM3/  CAPAl , CAPA2 , CAPBl , CAPB2 
C 

C  THIS  SUBROUTINE  EVALUATES  THE  FUGACITY  FOR  THE  SOLUTE  IN  SOLUTION 

C 

TERM1  =  Z+CAPBMIX* (1 .DO-DSQRT (2 .DO)  ) 

TERM2  =  Z+CAPBMIX* (1 .DO  +  DSQRT (2 .DO)  ) 

TERM3= (2* (YI*AI2+Y2*A2) /AMIX) - (B2/BMIX) 

TERM4=CAPAMIX/ (2 .D0*DSQRT (2 .DO) *CAPBMIX) 

TERM5=TERM4*DLOG (TERM1/TERM2) 

TERM6= (B2/BMIX) * (Z-1 .DO) 

TERM7=-DLOG (Z-CAPBMIX) 

DLNPHI=TERM6+TERM7+TERM3*TERM5 
PHI2H=DEXP (DLNPHI) 

C  WRITE (6, 100)  TERMl, TERM2, TERM3, 

C  WRITE (6, 100)  TERM4, TERMS, TERM6 

C  WRITE (6, 100)  TERM7, DLNPHI, PHI2H 

C  100  FORMATC  '  ,  5  (2X,  D15 . 6)  ) 

RETURN 

END 

C 

C 

C 


SUBROUTINE  SOLUEST ( T , P, Y1 , Y2 , Z , PHI2H, Y2EST) 
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IMPLICIT  REAL*8 (A-H, 0-Z) 

COMMON  /PARAM/  AlC , BIC , A2C , B2C 

COMMON  /PARAMl/  TCI , PCI , OMEGAl , TC2 , PC2 , 0MEGA2 , DEL12 
COMMON  /PARAM2/  A1 , A2 , A12 , AMIX, B1 , B2 , BMIX, CAPAMIX, CAPBMIX 
COMMON  /PARAM3/  CAPAl , CAPA2 , CAPBI , CAPB2 
C 

C  THIS  SUBROUTINE  ESTIMATES  THE  SOLUBILITY  OF  THE  SOLUTE  IN  THE 

C  SUPERCRITICAL  PHASE.  VALUES  FOR  THE  VAPOR  PRESSURE  OF  THE  SOLUTE 

C  AS  A  FUNCTION  OF  TEMPERATURE  MUST  BE  AVAILABLE.  ALSO  NEEDED  IS 

C  THE  SATURATED  LIQUID  VOLUME  AT  THE  TEMPERATURE  OF  INTEREST.  BOTH 

C  VAPOR  PRESSURE  AND  SATURATED  LIQUID  VOLUME  ARE  USED  TO  EVALUATE 

C  THE  POYNTING  FACTOR. 

C 

C  EVALUATION  OF  VAPOR  PRESSURE 

C 

RG=83 . 14D0 
BVP=7857 . 52D0 

VPLOG10= (-BVP/TC2) * (1- (T/TC2) ) / (T/TC2) 

VP=10.D0** (VPLOGIO) 

P2SAT=VP*PC2 

C  WRITE (6, 200)  P2SAT 

C  200  FORMAT ( '  ',2X, 'P2SAT  =',D15.6) 

V2SAT=438 . 18d0/2 . OdO 
C  WRITE (6, 201)  V2SAT 

C  201  FORMAT ( '  ',2X, 'V2SAT  =',D15.6) 
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Y2EST= (P2SAT/P) * (1 .D0/PHI2H) *DEXP (V2SAT* (P-P2SAT) / (RG*T) ) 
C  WRITE(6,202)  Y2EST 

C  202  FORMAT ( '  ',2X, 'Y2EST  =',D15.6) 

RETURN 

END 
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